Electron-positron interaction in light elements 
represented by atoms embedded in an electron gas 
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Mijnarends et al. [J. Phys. Condens. Matter 10, 10383 (1998)] contested the best existing 
calculations of positron annihilation rates in jellium and crystal lattices, pointing in this way at 
deficiencies of existing theories of electron-positron interaction in these materials. In the present 
work the local enhancement factors due to e"*" — e~ interaction in Li, Be, B, C, N and O are 
computed in a consequent many-body approach for core and as concerns lithium also for conduction 
electrons and compared to the results of existing approximations to this problem which avoid direct 
many-body calculations in metals, i.e. the local density, generalized gradient and weighted density 
approximations, as well as to experimental data. Conclusions about positron lifetime and — e~ 
correlation energy are also presented. Suggestions concerning annihilation rates in an electron gas 
agree with those of Mijnarends et al. 

PACS numbers: 71.60.-l-z, 78.70.Bj, 71.10.-l-x, 71.25.Cx 



I. INTRODUCTION 

The purpose of this work is a study of electron-positron 
interaction (EPI) in solids represented by an atom em- 
bedded in an electron gas. There are good reasons to 
believe that EPI in such a model gives valuable informa- 
tion on — e~ interaction. Indeed, in most calculations 
the local annihilation rates are assumed to depend only 
on the local electron density, in generalized gradient ap- 
proximation also on its local gradienti. In our work they 
are functionals of the whole electron distribution in an 
atom while the neighboring atoms are simulated by an 
averaged electron density and a constant positive back- 
ground. 

This model was used to compute the effect of EPI 
for light elements. However, the approach to electronic 
structure used in this work (i.e. presenting the pop- 
ulation of valence electrons by means of a single den- 
sity amplitude in the form proposed by Gondzik and 
Stachowial^) seems to hold only for alkalis^. So more 
or less reliable calculations of the enhancement of va- 
lence electrons in these elements were performed only for 
lithium. As concerns other elements, we had to limit our- 
selves to core electrons for which the density amplitude 
is equivalent to the wave function. 

The unknown effect of EPI on annihilation characteris- 
tics has always been considered as an obstacle in interpre- 
tation of annihilation data^. Nevertheless positron anni- 
hilation found many applications in studies of the solid 
state and in other domains (including biology, medicine 
and even marketing) and proved to give valuable infor- 
mation on electronic and ionic structuresS:-^. 

EPI in many electron systems has also been the subject 
of many investigations both theoretical and experimen- 
tal. A review of different approaches to this problem is 
given in Ref . H. 

In particular, the unknown effect of the interaction of 



the positron with core electrons (IPC) introduces an un- 
certainty which is an obstacle in the interpretation of 
experimental data in terms of properties of valence elec- 
trons. 

From the many attempts to estimate the effect of IPC 
(and EPI in general) on the annihilation data let us men- 
tion the local density approximation (LDA)2ii£, the gen- 
eralized gradient approximation (GGA)i, the weighted 
density approximation (WDA)ii and othersiSii^. 

In the worka^iii no attempt is performed to study the 
behavior of electrons in solids in presence of the positron 
basing on the equations of quantum mechanics. They all 
benefit of the results of jellium calculations. GGA takes 
into account the inhomogeneity of the medium, but has 
to introduce a phenomenological constant in order to get 
agreement with experimental measurements of positron 
lifetimes in metals. This agreement was contested as con- 
cerns other applications by Mijnarends et al.^'* who per- 
formed two-detectors Doppler broadening measurements 
in Al. WDA has the merit of enforcing a total charge 
of the screening cloud around the positron equal to one 
electronic charge. 

Among the many-body calculations the most advanced 
are those of Sormann^^. But even in them no attempt to 
reach self-consistency is undertaken. 

In the present work we will develop the approach pro- 
posed by Gondzik and Stachowiak^ in order to treat a 
positron in an electron gas and called mnemotechnically 
HNC (hypernetted-chain) following in that the work of 
Kallio et alpi^. The approach of Ref. is particularly 
simple, and leads to reasonable results. So one can hope 
that generalizing it to real solids will be easier than in 
the case of other approaches. The basic equations used 
in the present work have been proposed in Ref. ^3 Etnd 
the method of solving them for the anisotropic case (i.e. 
for the positron beyond the center of symmetry) is pre- 
sented in Ref. 0. This work constitutes an application 
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of the methods elaborated in these last papers. 

Indeed, some information about e"*" — e~ interaction in 
lithium was already obtained along this lineiS*ii. How- 
ever, in these works we were unable to solve the occur- 
ing integro-diflerential equations in the two dimensions 
needed if the spherical symmetry is broken (it will be bro- 
ken for the positron outside the nucleus). In the mean- 
time a method to solve such equations was elaborated 
and then applied to Iithiumi24i2i2£. 

In this work we had to limit ourselves to light ele- 
ments, since at present we are able to deal only with 
core electrons in the s orbital state. So we will perform 
calculations only for Li, Be, B, C, N and O. 

The present state of the art can be characterized best 
basing on the recent enlightening work of Mijnarends 
et ali^. These authors interpret two-detectors Doppler 
broadening studies of Al in the following way. 

They assume three models of e+ — e~ interaction which 
they label LDA, LDA' and GGA. The first two models 
apply the local density approximation to e+ — e~ inter- 
action using the results of electron gas theory, in the first 
case those of Arponen and Pajanne (AP)^, in the second 
case those of Lanttc^^. The third case corresponds to cor- 
recting the results of the local density approximation by 
adding a gradient correction to the enhancement factors 
of AP according to the work of Barbiellini et al.A. Mi- 
jnarends et al. find that the LDA' model gives the best 
agreement with experiment, though a small correction 
by applying a lesser than in Ref. 1 gradient correction 
could maybe lead to some improvement. These results 
illustrate the present state of knowledge of — e~ inter- 
action both in an electron gas and in real metals. Com- 
parison of GGA with experiment performed in Ref. |3 
led to disagreement. 

As concerns deviations of enhancement factors from 



the local density approximation in metal lattices, this 
problem according to Ref. Q is of lesser importance. But 
this conclusion follows from experimental considerations 
and means that a satisfactory theory of e"*" — e~ inter- 
action in metal lattices is not existing in spite of more 
than thirty years of research in this direction. Moreover, 
the LDA' model is based on the calculations of Lantto 
which assume an oversimplified trial function of the Jas- 
trow type. This trial function neglects as well momen- 
tum dependence of e+ — e~ scattering and dependence of 
— e~ correlations on the distance from the positron - 
effects well established in physics and included in other 
calculations2iiSi2^iS^. 



II. COMPUTATIONS 

The model 

Applying the approach of Gondzik and Stachowiak^ to 
the problem of positron screening in an electron gas, we 
describe the electronic structure of an atom with two elec- 
trons core embedded in jellium with two functions ijji{r) 
and ^2{i')- "tpiif) is the wave function of core electrons 
and '02(?') is the density amplitude of valence electrons. 
■!/'2(^) is equal to the density of conduction electrons in 
the model. These two functions obey the appropriate 
Kohn-Shamlike equations (in atomic Hartree units which 
will be used throughout the paper): 



[-iV2-f F(r)]V2(r)= 



(1) 



where 



J 



dr' 



, ^?(r') 



>i(r')-d(r') 



+ VHL{2^f{r) + ^Jl{r)} - VhlIpo}- 



(2) 
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Z is the charge of the nucleus, d(r) is the distribution of 
the positive charge in the electron gas. It is equal 



d(r) 



Po 




for 
for 



r > Rws, 
r < Rws 



where 



3(Z-2) 
ATT{Rwsf 



(3) 



(4) 



Rws is the radius of the Wigner-Seitz sphere. Eqs. ^ 
- ^ limit us to elements having a two-electrons core. 
Vhl{p} is the Hedin-Lundqvist exchange-correlation cor- 
rection for an electron gas of density p^. The Lagrange 



multiplier E2 which should occur on the right-hand side 
of the second equation ^ is normalized to zero by the 
last term in the formula ((2Jl , while Ei is the energy eigen- 
state of core electrons. 

We admit, of course, that from the elements treated in 
this work only lithium can be described as above^iii. But 
since core electrons are only weakly aff'ected by valence 
electrons, the way of presenting the electronic structure 
of these lasts is of little importance for our purpose. In 
fact in this work we used different ways of approximating 
the electronic structure of valence electrons, but this had 
a negligible infiuence on IPC. 

D was chosen in order to satisfy the obvious normal- 



3 
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FIG. 1: Comparison of the density of valence electrons for a 
lithium atom in an electron gas - full curve for D — 1.24 and 
dotted curve for D = 1 and the analogous density (averaged 
over direction) for metallic lithium in HNC approximation- 
dashed curve. 



Hum where the value D = 1.1 was obtained. In Fig. 
the distribution of valence electrons V'K'') is shown for 
D = 1.24 and compared to the distribution of valence 
electrons in metallic lithium obtained in HNC formalism^ 
after averaging over directions. It is visible that the value 
D — 1.24 reproduces better the electron distribution in 
the immediate neighborhood of the atom. 

Introduction of the positron 

We benefit of the result obtained according to the the- 
ory of liquids by Kallio et al.^^ and concerning the role 
of the light mass of the positron in EPI. This suggests to 
describe the electronic structure of the model in presence 
of a positron at r^ by means of the equations 



Xi(rp,r) = Ei{rp)xi{rp,r), 



lv^ + V{r) + ^Wirp,v) 



X2(rp,r)=0 (6) 



ization condition 
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Jo 



Z -2. 



(5) 



This requirement could be satisfied for lithium where the 
value D ~ 1.24 leads to satisfying Eq. (0 and in beryl- 



where the screened electron-positron potential W{rp,r) 
is defined as 



W{rp,r) 
Here 



1 



r - r„ 



Wpirp,r) + W^c{rp,r). (7) 



M^p(rp,r) =2 / dr' 



, x?(rp,r')-^?(r') 



dr' 



, xi(rp,r')-Vi(r') 



(8) 



WU^^p, r) = yHL{2x?(rp, r) + xli^p. r)} - VuLi'^^Hr) + ^|(r)}. 



(9) 



r 



Xi indicate the form of the functions ijji in presence of 
the positron. The positron distribution is described, of 
course, by the positron wave function. 

Since the dependence of IPC on the perturbation of va- 
lence electrons by the positron is negligible, the equation 
for xi can be solved separately. However, we consider as 
a better approximation to solve exactly the equations © 
for the easy problem of the positron on the nucleus and 
to compute xi from the equation 



l^^ + V\v) + \w\vp,v) 



Xi(rp,r) 
= E^{rp)xi{vp,v). (10) 



Since the positron in metals is always screened by valence 
electrons, we found that when computing xi it would be 
appropriate to freeze the valence electrons in the state 



they acquire for the positron on the nucleus. This will 
allow to determine the potentials and W"^ in Eq. H1U|I . 
This equation will be solved afterwards for all values of 

Yp. 

Let us compute the solution of Eq. (0 for the positron 
on the nucleus, assuming that the core is frozen. We get 
in this way the potential 

g(r)= ('-O + VHL{2i^l{r) + xi(0,r)} 



'VHL{2i^l{r)+i4{r)}- 



(11) 



The potentials in Eq. pO(l can now be written in the 
form 



V\r) = V{r) + ig(r), 



(12) 
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r - r„ 



dr' 



,X?(rp,r')-V?(r') 



(13) 

Solution of the equations 

From test calculations it follows that the effect of 
positron interaction with conduction electrons depends 
only slightly also on the polarization of the core by the 
positron. So the equation for X2 can be solved separately. 

The functions Xi ^'^^ presented in the form 



X^(rp,Y) = A,{rp)e + r,(rp, r) 
where the function 



(14) 




FIG. 2: The density amplitude of valence electrons in lithium 
as seen by positrons x^iTpyrp) (full curve) compared to the 
true density amplitude ^/'2(rp) (dashed curve). 



n(i-p,r) 



n=0 



^n{rp,r)Pn{cOSd) 



(15) 



is devoid of the cusp occuring in Xi ^-t the positron. 
P„ (cos 1?) are Legendre polynomials where d is the an- 
gle between r and rp, s =| r — rp |. From the assumption 







(16) 



(the simplest one but not necessarily the only possible) 
it follows that a — 1/2. 

We have to limit ourselves for technical reasons to two 
terms in the expansion H15|l . This leads to some problems 
especially for higher values of rp. They are discussed in 
more detail in Ref. 18. 

The contact densities of electrons on the positron are 
provided by the Ai{rp) coefficients. An important infor- 
mation is also contained in the energy eigenvalue Ei{rp). 
Since 2Ei{rp) will enter the positron Hamiltonian as a 
contribution to the positron potential, this quantity can 
be interpreted as effective attraction (of chemical char- 
acter) between the positron and the nucleus due to col- 
lectivization of core electrons. 



III. RESULTS 

We were able to obtain numerical values for local an- 
nihilation rates of positrons in lithium. The results are 
shown on the figures. For other light elements we were 
able to study the effect of the interaction of the positron 
with atomic cores. 

We call LDA the local density approximation, men- 
tioning each time what results are used for describing 
the properties of a homogeneous electron gas: the ones 
of Gondzik and Stachowiak^ or PHNC^'^ (eventually the 
formula of Bororiski and Nieminen (BN)^) . Since the in- 
teractions of the positron with conduction and with core 
electrons can be to a large degree considered as indepen- 
dent, we introduce the term local partial density approx- 
imation (LPDA) which means that only core electrons 



or only conduction electrons (depending which electronic 
state is studied) contribute to the local density. 

Note that the enhancement factors following from Eqs. 
© are approximate as concerns conduction electrons 
since they neglect the difference of scattering on the 
positron (and also on the nucleus) by different elec- 
tronic states. This is why it is appropriate to compare 
them while using LPDA with the results of the Gondzik- 
Stachowiak approach. This problem does not occur in 
the case of core electrons. 

The enhancement factors for core {i = 1) and for va- 
lence electrons {i — 2) are given by the square of the 
enhancement amplitude Wi which is defined as 



'w^(rp) 



(17) 



On Fig. 12 the density amplitude as seen by positrons 
X2{fp,rp) is compared to the real density amplitude 
i'2irp). 

The enhancement amplitude of conduction electrons 
calculated according to IjlTI) agrees quite well with LPDA 
predictions at rp = 0. Also for rp much bigger than the 
position of the node in tp2 it approaches the LPDA value. 
Note, however, that for big r's its value is higher than ex- 
pected on ground of the local density approximation. We 
can interpret this effect as a result of increasing electron 
density between (approximately) the position of the node 
in 'ip2 and r < 2.5. This result is confirmed by the recent 
calculations of Boronski and Stachowiak^ who study the 
possibility of applying a grid method for problems of that 
kind. A similar effect has been observed for vacanciesi^ 
where a maximum of the enhancement factor occurs in 
the region of increasing electron density. In the present 
case we rather attribute this effect to the lack of period- 
icity of the model. In the intermediate region A2{rp) is 
greatly affected by the displacement of the node (from 
0.9 to 0.76 a.u.) due to interaction with the positron. We 
consider this effect as very important, impossible to ob- 
tain using other approaches to EPI. Such an effect has 
been observed experimentally by Chiba in MgO^. 
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The enhancement amphtude for core electrons is pre- 
sented in Fig. for Li. Remark that in the immediate 
vicinity of the nucleus the enhancement amplitude is a 
little bigger than LPDA predictions, while for higher val- 
ues of r it falls below the LPDA curves. 

Fig. Oh shows that the effective enhancement for 
lithium defined as 



"eff 



(rp) 



(18) 



reproduces quite well LDA predictions. A2{rp) has been 
obtained from A2{rp) by renormalizing it, taking account 
of the increasing error while increases when solving 
Eqs. © and the deviation of the approximation of Ref.Q 
from the exact value assumed to correspond to PHNC— . 

Of course, the local enhancement factor presented in 
Fig. |3Ji as well as the local enhancement of conduction 
electrons influence not only the positron lifetime but also 
angular correlation results. 

As concerns the energy of e"*" — e~ correlation , one 

can assume that conduction electrons give to it a constant 
(independent of r) contribution Ey, while the contribu- 
tion of core electrons is provided by 2Ei{rp). In Fig. 0^ 
we compare = 2Ei{rp) + Ey to the appropriate pre- 
dictions of LDA for the energy of e+ — e~ correlation. 
The jellium values were chosen according to PHNC^^. 

In Fig. we compare the positron wave function in 
lithium calculated using the LDA approximation for the 
e"*" — e~ correlation potential to the wave function ob- 
tained when the effect of correlation was estimated ac- 
cording to our results. 

We found also total annihilation rates for these two 
cases. We used the LDA approximation for the enhance- 
ment factor estimating that in the light of Fig. [Jb it is 
quite justified. In the first case, when we applied the LDA 
approximation to the PHNC enhancement factor (using 
the PHNC formulaiS^) and to the e+ — e~ correlation po- 
tential (with the formula from Ref. I29t) we obtained a 
total annihilation rate of 3.63 10^/s. The electron densi- 
ties were calculated on base of the FLAPW (full poten- 
tial linearized plane wave) method (the corresponding 
numerical code was WIEN9£'^"). 

The bigger penetration into the core region by the 
positron as shown in Fig. 0j3 results in slightly higher 
annihilation rates than in LDA because of high electron 
densities inside the core. Thus, for the potential calcu- 
lated according to this work, the corresponding rate was 
3.68 lOVs. 

The calculations performed with the BN 
enhancement^^ and the LDA correlation potential 
give a rate of 3.36 10^/s. On the other hand, the 
calculations performed with the BN enhancement and 
our correlation potential give a rate of 3.42 10^/s, what 
is accidentally in perfect agreement with experiment. 
The experimental annihilation rate for lithium is equal 
3.436 107a22. 

When calculating the positron wavefunction we used 
the numerical program by M. Puska'^^ based, in princi- 



ple, on superposition of atomic densities. We changed it 
slightly, however, in order to have the possibility to use 
electron densities and potentials calculated within the 
FLAPW method. Moreover, we generalized it on an ar- 
bitrary non-LDA e+—e~ correlation potential. Note that 
Puska himself obtained an annihilation rate of 3.28 
basing on electron densities calculated according to the 
LMTO-ASA (linearized muffin-tin orbital - atomic sphere 
approximation) approach. Sormann and Sob^ remarked 
already that different ways of computing band structure 
lead to predicting different positron annihilation charac- 
teristics. 



IV. CONCLUSIONS 

In this work we tried on the simple example of lithium 
to find in a consequent many-body way (by solving in a 
non spherically symmetrical surrounding the appropriate 
nonlinear integro-differential equation derived on ground 
of the theory of liquids) deviations of local positron an- 
nihilation rates in a metal lattice from the local den- 
sity approximation. We came to the conclusion that in 
lithium no drastic deviation from LDA occurs as con- 
cerns the annihilation rate. However, as concerns valence 
electrons, we observe a displacement of the node from 
0.9 a.u. to 0.76 a.u., similar to that detected by Chiba 
in MgO^. We found also a deviation of the energy of 
electron-positron correlation from the local density ap- 
proximation usually assumed in calculations. We found 
instead a term in the energy which we interpret as at- 
traction of chemical character between the nucleus and 
the positron due to collectivization of core electrons. 

So we computed the total positron annihilation rate 
in lithium using the local density approximation for the 
local annihilation rates. We used the numerical ATSUP 
code of Puska™ for calculating positron wave functions, 
however, unlike Puska, we performed band structure cal- 
culations using the FLAPW code labelled WIENgS^o. 
Moreover we added to the positron potential the term 
describing the positron-nucleus attraction as obtained in 
our calculations. As concerns jellium annihilation rates. 



we used the formula of Borohski and Niemine: 



In this 



way we obtained perfect agreement with experimental 
values and with the conclusions of Mijnarends et ali^ 
concerning the LDA' model. 

We confess that we are not very happy about this re- 
sult. The formula of Borohski and Nieminen is the most 
frequently used in calculations trying to reproduce ex- 
perimental data. Obviously it was found to give the best 
agreement. However, it is based on the calculations of 
Lantto which assume a Jastrow type trial function. This 
function neglects the momentum dependence of electron- 
positron scattering and the influence of the positron on 
electron-electron correlations. So we would have more 
confidence in calculations which take these effects, so well 
established in physics, into account. However, the anni- 
hilation rates obtained from them are obviously too high 



6 



3.0n . . . . 1 12 




0.o4 ■ 1 ' 1 ' 1 "~ i I 0^ ■ 1 ^ 1 ^ 1 ■ 

0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

rp(a.u) rp(a.u.) 

FIG. 3: a) The enhancement amphtude of core electrons in Li according to this work (full curve). The dashed curve shows the 
LPDA prediction. Additionally, the core electron distribution A-Kr'^^i(r) is plotted in arbitrary units (dashed-dotted curve), 
b) Effective enhancement e(r-p) inside an atom of lithium embedded in an electron gas calculated in this work (solid curve) and 
according to LDA (dashed curve), GGA (dotted curve) and WDA (dashed-dotted curve). 




FIG. 4: a) Comparison of (r^) — 2Ei(rp) + Ev in Li 

(full curve) and the energy of electron-positron correlation 
according to LDA (dashed curve). The energies have been 
normalized in such a way as to coincide at Vp = 2.1 a.u. 
b) The positron density in Li along the (100) direction. The 
dashed curve was obtained while using the e"*" — e~ correla- 
tion energy in LDA approximation. The full curve follows 
from including the effect of nucleus-positron attraction due 
to collectivization of core electrons. 
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FIG. 5: Comparison of different formulas for the positron 
annihilation rate A(rs) in an electron gas (dashed curve 
- Gondzik-Stachowiak", dashed-dotted curve - Barbiellini 
et a^, solid curve - PHNCS^, dotted curve - Bororiski- 
Nieminen«) . 



(Fig. O). So we have to accept the Borohski-Nieminen 
formula as an expression poorly explained theoretically, 
but describing pretty well experimental data. 

Among other results contained in the paper let us men- 
tion calculations of positron-nucleus attraction for Be, 
B, C, O and N due to collectivization of core electrons, 
of core enhancement for these elements (Appendix 



and description of the distribution of the electronic cloud 
screening the positron in different positions (Appendix 
0. 
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APPENDIX A: POSITRON-CORE 
INTERACTION IN LIGHT ELEMENTS 

For light elements heavier than beryllium the approach 
described by Eqs. and (0) failed. Indeed, in boron it 
yielded D = 0.59 and collapsed in carbon where the con- 
dition lO needed a value of D equal 2.5 • 10"'^. This, of 
course, is connected with the properties of valence elec- 
trons which are no longer itinerant. For this reason the 
condition Q was replaced for boron, carbon, nitrogen 
and oxygen by the condition 



in 



cRws 



(Al) 

while in Eq. lAip the role of Rws is played now by 
cRws- The normalization condition for is the usual 
one. Note that the dependence of IPC characteristics on 
the values of D and c in Eqs. 0] and ljAl|l is negligible. 
Our calculations show that the enhancement amplitude 
for core electrons in Be, B, C, O and N is even lower 
from the LDA curve than the corresponding figure for 
Li. In general, however, its behavior is quite similar to 
that shown in Fig. Unfortunately, the effective en- 
hancement for those elements could not be computed for 
reasons explained in Section II. 

As concerns the energies of positron-electron correla- 
tion for Be, B, C, O and N the effect of core electron col- 
lectivization leads to an attraction between the positron 
and the nucleus. The effect is bigger for increasing atomic 
mass of the element. 



1.8 . Fig. Ul shows the behavior of the screening cloud 
X2(rp,r) — "02 (i") foi' t'^s same values of r and r^. One 
should remember that the corresponding density ampli- 
tude changes sign when crossing the node, and this fea- 
ture persists in the presence of the positron. Moreover, 
the node is shifted by the positron. It is striking that the 
highest density of the screening cloud occurs close to the 
nucleus even for the positron well beyond the node region 
of the density amplitude. This suggests that our figures 
are less reliable for the positron in the interstitial region, 
since the screening cloud extends over several atomic 
cores and these facts are not included in our model. 




s (a.u.) 



FIG. 6; X2(i'p, s -I- Vp) — 'ip2{s + rp) in Li for rp equal a) 0.3 
b) 1.8 . s is taken along the line connecting the nucleus and 
the positron. 




APPENDIX B: THE DENSITY AMPLITUDE 

X2(rp,r) 

Fig. shows the change of the density amplitude 
X2(rp,r) — ■02(1") for r along the line connecting the 
nucleus and the positron and rp equal a) 0.3 and b) 



FIG. 7: The screening cloud distribution in Li around a 
positron at rp equal a) 0.3, b) 1.8 . s = r — is taken along 
the line connecting the nucleus and the positron. 
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